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The study of a certain class of matrix integrals can be motivated by their interpretation 



^ ■ as counting objects of knot theory such as alternating prime links, tangles or knots. The 

r-| \ simplest such model is studied in detail and allows to rederive recent results of Sund- 

I . berg and Thistlethwaite. The second non-trivial example turns out to be essentially the 
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1. Introduction 

Using random matrices to count combinatorial objects is not a new idea. It stems 
from the pioneering work of [0 , which showed how the perturbative expansion of a simple 
non-gaussian matrix integral led, using standard Feynmann diagram techniques, to the 
counting of discretized surfaces. It has resulted in many applications: from the physical 
side, it allowed to define a discretized version of 2D quantum gravity and to study 
various statistical models on random lattices IQJ^]. From the mathematical side, let us cite 
the Kontsevitch integral P,^J^] , and the counting of meanders and foldings PJ|] . 

Here we shall try to apply this idea to the field of knot theory. Our basic aim will be 
to count knots or related objects. In the next sections, we shall define these objects; then 
will follow a brief overview of matrix models and how they can be related to knots. We 
shall then explain the counting of alternating links (following |]T^) in section 4; study a 
generalized model (the ABAB model [|ll|) in section 5, which will lead us to digress and 
consider summations over Young tableaux; and finally, discuss further generalizations in 
section 6. 



2. Knots, links and tangles 

Let us recall basic definitions of knot theory (from a physicist's point of view; the 
reader is referred to the literature for more precise definitions) . A knot is a smooth circle 
embedded in IR . A link is a collection of intertwined knots. Both kinds of objects are 
considered up to homeomorphisms of IR . Roughly speaking, a tangle is a knotted structure 
from which four strings emerge. 

In the 19*^^ century, Tait introduced the idea to represent such objects by their pro- 
jection on the plane, with under/over-crossings at each double point and with minimal 
number of such crossings (fig. 1). 




6i 



Fig. 1: A reduced diagram with 6 crossings. 



We shall now consider such reduced diagrams. To a given knot, there corresponds a 
finite number of (but not necessarily just one) reduced diagrams. We shall come back later 
to the problem of different reduced diagrams which correspond to the same knot (or link, 
or tangle). 

To avoid redundancies, we can concentrate on prime links and tangles, whose diagrams 
cannot be decomposed as a connected sum of components (Fig. 2). 




Fig. 2: A non-prime diagram. 



A diagram is called alternating if one meets alternatively under- and over-crossings as 
one travels along each loop.u From now on, we shall concentrate on alternating diagrams 
only, since they are easier to count. There are two reasons for that. 

The first reason is that there is a simpler way to characterize if two reduced alternating 
diagrams correspond to the same knot or link (than the general Reidemeister theorem ||12|| ). 
Indeed, a major result conjectured by Tait and proven by Menasco and Thistlethwaite [^ 
is that two alternating reduced knot or link diagrams represent the same object if and only 
if they are related by a sequence of moves acting on tangles called "fiypes" (see Fig. 3). 





Fig. 3: The flype of a tangle. 



The second reason is that there is a correspondance between alternating diagrams and 
planar diagrams (see for example flj]), which will be explained now as we discuss matrix 
integrals. 



Even though it may not seem obvious, there are knots that cannot be drawn in an alternating 
way - starting with eight crossings. 



3. Matrix integrals 

Let us now start from a completely different angle and consider the following matrix 
integral: 

Z^^\g)= / dMe ^ ^ ^ ^ (3.1) 



M is a A^ X A^ hermitean matrix; (7 is a real parameter, which should be chosen negative 
to make the integral convergent. 

As an application of Wick's theorem, the perturbative expansion of Z^^' in powers 
of g can be made using the following Feynman rules: one should count all diagrams made 
out of vertices "^^ = gN5qi5jk and propagators {MijMj^^Q = )= =. [ = jj5u5jk. Due 

1 * 

to the double lines, these diagrams form so-called fat graphs which can be identified with 
triangulated surfaces. Each diagram has a weight 



symmetry factor 



where V, E, F are the number of vertices, edges, faces of the triangulated surface (the 
factor A^^ comes from the summation over internal indices). The symmetry factor (i.e. 
the order of the automorphism group of the diagram) is of little importance to us and 
we shall not discuss it any further. Note that the power of A^ is simply A^^ where x is 
the Euler-Poincare characteristic of the triangulated surface. If we take the logarithm, 
which amounts to considering only connected surfaces, then we have the following genus 

expansion: 

00 

where Fh is the sum over surfaces of genus h. In particular, if we consider the large A^ 
limit, we see that 



^1 



N^oo A -^-"^ symmetry factor 

planar graphs 

is the sum over connected "planar" diagrams (i.e. with spherical topology). F{g) is the 
quantity we are interested in. The formal power series F{g) = ^ fpgP turns out to have, 
as is well-known, a finite radius of convergence (which allows to analytically continue it 
to positive values of g, as will be explained later). The position and nature of the closest 



singularity Qc determines the asymptotics of /p as p ^ oo i.e. of the number of planar 
diagrams with large numbers of vertices. 

In order to connect with knot theory, we take any planar diagram and do the follow- 
ing: starting from an arbitrary crossing, we decide it is a crossing of two strings (again 
there is an arbitrary choice of which is under/over-crossing) . Once the first choice is made, 
we simply follow the strings and form alternating sequences of under- and over- crossings. 
The remarkable fact is that this can be done consistently (Fig. 4). If we identify two alter- 
nating diagrams obtained from one another by inverting undercrossings and overcrossings, 
then there is a one-to-one correspondence between planar diagrams and alternating link 
diagrams. So the function F{g) also counts alternating link diagrams with a given number 
of crossings. 





Fig. 4: A planar diagram and the corresponding alternating link diagram. 



A more detailed discussion of the properties of the resulting link diagrams will be made 
in the next section. For now, we shall address the question of the number of connected 
components of the link (as a 3-dimension object). Indeed, there is no reason for the diagram 
to represent a simple knot, and not several intertwined knots. In order to distinguish them, 
we introduce a more general model, which we shall call the intersecting loops 0{n) model. 
If 71 is a positive integer, then consider the following multi-matrix integral: 



Z^''\n,g)= /HdM^e 



A^tr (-i ELi Ml + f E;6=i MMbMaM, 



and the corresponding free energy 

F{n,g) 



lim 



log Z(^\n,g) 



(3.2) 



(3.3) 



This model has an 0(n)-invariance where the Ma behave as a vector under 0{n). Its 
Feynman rules are a bit more complicated since we should draw the diagrams with n 
different colors. The colors "cross" each other at vertices just like strings in links (Fig. 5). 




Fig. 5: A planar diagram with 2 colors. 

So what we have done is allow each loop in the link to have n different colors. This 
is in itself an interesting generalization of the original counting problem. Indeed, we can 
write: 

oo 

F{n,g) = Y,F''{g)n'^ (3.4) 

fc=i 

where F^{g) is the sum over alternating link diagrams with exactly k intertwined knots. 
We see that the links are weighted differently according to their number of connected 
components. 

But there is more. The expression (|3.4| ) is an expansion of -F(n, g) as a function of n 
around 0; it provides a definition of F{n, g) for non-integer values of n. In particular, we 
have the following formal expression 

^, dFM 

on |?i = 

for the sum over alternating knots (this is the classical replica trick). Therefore, if one 
computed F{n, g) for arbitrary (non-integer) values of n, one would have access to the 
generating function of the number of alternating knots. Of course, it might seem difficult 
to solve our model for all n; we shall discuss this again in the conclusion. 



4. The one matrix model and the counting of links 

Let us now come back to the one-matrix model and show how one can derive explicit 
formulae for the counting of prime alternating links. We recall the partition function 



iM^ I iVtr(-iM2 + fM4) 

Z^^\g) = / dMe ^ ^ ^ ^ (4.1) 



and the corresponding free energy 



F(^)= hm i^logZW(r7) (4.2) 



We also define the correlation functions: 

G2n{g) = lim (^trM^-) (4.3) 

Whereas the perturbative expansion of F{g) generates closed digrams (and therefore al- 
ternating links), the G2n(fi') count diagrams with 2n external legs. In particular, we shall 
be interested later in r{g) = G4^{g) — 2G2{g)'^ which counts connected diagrams with 4 
legs, i.e. alternating tangles. 

There are various methods to compute all these quantities. We shall briefly recall the 
simplest one: the saddle point method. 

4-1. Saddle point method for the one-matrix model 

We start from Eq. ( |4.1| ) and notice that the action and measure are C7(A^)-invariant; 
therefore we can go over to the eigenvalues A^ of M: 

Zi9) = / n dA.A[A.]2 e'' ^'^^ ^ ^^^ ^ ^^^ ^ (4.4) 

up to an overall constant factor. Here A[-] is the Van der Monde determinant. In Eq. ( [4.4| ), 
the action is of order A^^ while there are only A^ variables of integration. Therefore, in 
the large A^ limit, a saddle point analysis applies. It is easy to see that as A^ — i> oo, the 
eigenvalues A^ form a continous saddle point density p(A) whose support is an interval 
[—2a, 2a]. In order to solve the saddle point equations, it is convenient to introduce the 
resolvent 



c<;(A) = lim ( — tr 



N^oo \N X- M 

" dA' "^'"> 



(4.5) 



-2a A — A' 

Then the saddle point equations read: 

uj{X + iO)+u{X~iO)-X + gX^ = A G [-2a, 2a] (4.6) 

This is a simple Riemann-Hilbert problem which can be solved: 

uiX) = h- ^gX' - (^-^gX' + ^ - ga'^ ^ X^ - Aa^ (4.7) 



with a^ = I — — -. What we have found is the generahzed semi-circle law (for g 







we recover the usual GUE). Since a; (A) is a generating function of the G2n, we can extract 

1 



Also, we find 



G2{g) 



V{g)=G^{g)-2G-,{g) 



1 




§a2(4-a2) 



a^(a2-l)(2a2-5) 



F{9) = - log a 



1 



(a2-l)(9-a^) 



2 ° 24 

Note that all these expressions can now be analytically continued to (7 > all the 
way to the singularity Qc = 1/12. This has a simple interpretation: changing the sign of 
g > corresponds to making the potential in which the eigenvalues lie unstable; however, 
there is still a local minimum at the origin and since the large A^ limit is a classical limit, 
the eigenvalues cannot quantum tunnel to the unstable region and therefore remain in the 
valley (Fig. 6). However, as g reaches its critical value gc, the eigenvalues begin overflowing, 
which causes the singularity. 




Fig. 6: The potential for g > 0; analytic continuation is possible as long as tiie 
eigenvalues stay trapped inside the central valley. 

Of course, F{g) is not yet the counting function of prime alternating links. There are 
two separate problems to resolve: 

1) Are the diagrams reduced (i.e. do they have minimal crossing number)? Do they 
correspond to prime links? 

2) What about the flype equivalence? One should count only once different diagrams 
which are flype-equivalent. 

We shall address them now. 



4.2. Primality and Minimality 

The diagrams obtained from the matrix model can have "nugatory" crossings or "non- 
prime" parts (Fig. 7). 



D_ . _C^ . _ 



Fig. 7: First terms in the perturbative expansion of the 2-point function. 

Note however that aU these unwanted features appear as part of the two-point function. 
Therefore, in order to remove them, we must simply set the two-point function equal to 1! 
This is achieved by introducing an additional parameter t in the action: 

,^,, r iVtr(-|M2 + fM4) 

Z^^\t,g)= / dMe ^ ^ ^ ' (4.8) 



Of course, t can be absorbed in a rescaling of M, so the model is essentially unchanged. 
However we can now ask that t be chosen as a function of g such that 

G2itig),g) = l (4.9) 

We can solve this equation; the auxiliary function a{g) introduced earlier is now the solution 
of a third degree equation 

27g={a'^ -l){4-aY (4-10) 

equal to 1 when g = 0; and t{g) is given by 

tig) = la'{g)iA-a'{g)) (4.11) 

The function r{g) := r(t{g),g) is then the counting function for reduced alternating tan- 
gle diagrams. Similarly, F{g) defined by ■^F{g) = jG4{t{g),g) counts alternating link 
diagrams. We find in particular that the singularity of F{g) = ^ fpg^ (given by the 
equation gc/t^{gc) = 1/12) has moved to gc = 4/27; taking into account the power of the 
singularity, we find that the rate of growth of the number of alternating diagrams with p 
crossings is 

fp ^~°° const 6.75Pp-5 (4.12) 



A similar result was found in |T^ 



4-3. Flype equivalence 



The more serious problem we have to resolve is that we are not really counting links: 
we are counting diagrams, and links are flype-equivalence classes of diagrams. Here we 
shall follow 



Let us take a closer look at the action of a flype (Fig. 3). The key remark is that 
it acts on tangles (i.e. four-point functions), but more precisely on two-particle reducible 
(2PR) tangles. This leads naturally to the idea of introducing skeleton diagrams: a general 
connected tangle can be created by putting two-particle-irreducible (2PI) diagrams in the 
"slots" of a fully two-particle-reducible skeleton diagram. We then expect that the 2PR 
skeleton will be modified by the flype-equivalence, whereas the 2PI pieces (or more precisely 
the corresponding skeletons, see below) will be unaffected. Let T{g) = G4{g) — 2G2{gY 
be the counting function of connected tangles and D{g) of 2PI tangles; then T{D}, that 
is the power series obtained by composing V{g) and the inverse of D{g), is the counting 
function of fully 2PR skeleton diagrams (Fig. 8). 



\ 




\ 



/ 



D 
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Fig. 8: The whole set of diagrams F built out of the 2PI diagrams D. 



i-r(3) 



It easy to see from general combinatorial arguments that D{g) = T{g) 1,^(1 and 
therefore 



1 



T{D} = - ^1 - D - y/{l-D)^-4D 



(4.13) 



Inversely if D{g) = g -{- C(fi')7 then ([F] is the counting function of fully ^P/ skeletons 
diagrams (Fig. 9). 

From the solution of the one-matrix model, one obtains 



C[r] 



i + r 



2 11 

+ 2-r 



2(r + 2)3 



1 + lor - 2r2 - (1 - 4r) 



3/2 



(4.14) 



As we mentioned earlier, after taking into account the flyping equivalence, Eq. ( |4.13| ) 
will be modified, but not Eq. (^4.14| ). To see how it works, let us show how the counting of 
Fig. 8 is redone (Fig. 10). 




Fig. 9: The set of 2PI diagrams built out of general diagrams F (plus the single 
crossing g). 
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Fig. 10: Taking into account the flyping equivalence forces us to distinguish simple 
crossings from non-trivial 2PI diagrams (marked with a circle). There is only one 
term g(^ because the other term is obtained by a flype. 

More generally, we can redo the simple combinatorics to find the generating function 
of the 2PR skeletons, but this time taking into account the flype equivalence. We find 



ngX} 



1 



(i+g-c)~\ {i~9+cr-8c 



9' 



1-9 



This is to be combined with the (unaltered) matrix model data 



C[r] 



+ 2-r- -- 



1 



1 + lor - 2T^ - (1 - 4r) 



3/2 



1 + r ' ' 2(r + 2)3 

In practice, this means that T{g) is given by an implicit equation: 

t{9) = r{9,m9)]} 



(4.15) 



(4.16) 



(4.17) 



which can be reduced to a fifth degree equation. From the generating function of tan- 
gles T{g) we can go back to the generating function of closed diagrams F{g); we find in 
particular that the singularity has been displaced again, so that ii F{g) = Xl^o fp9^^ then 



jp ~ const 



101 + V2l00iy _7/2 
40 j ^ 



(4.18) 



where (101 + v/2T00l)/40 = 6.14793 . . . This result was first obtained in |T| . 
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5. The ABAB model and character expansion 

We shall now inspect the n = 2 case of the general 0{n) model (|3.2| ). There are various 
reasons that this model is of particular interest, and we shall discover some of them along 
the way. Let us rewrite the partition function 

r iVtr i-UA^ + S2) + f (A^ + B^) + ^{ABf) 



We see that we could introduce two coupling constants a and (3: 



tr (-1(^2 + S2) + f (A^ + B^) + f (AS)2) 



r ^ 

zfj^AB(a,l3) = jdAdBe ^ ' (5.2) 

(which amounts to introducing "interaction" between the two colors of loops). For a = 
(3 = g we recover the 0(2) model. The more general model with a and /3 arbitrary is not 
necessary for the original counting problem, but since it turns out that we can solve it 
equally easily, we shall keep the two coupling constants. Note that when a ^ (3 the 0(2) 
symmetry of the model is broken. This is even more apparent if we make the change of 
variables X = ^{A + iB), X^ = ^{A- iB): 

.^^ f , Artr(-XXt+6X2xt2 + |(XXt)2 + |(x4 + Xt4)) 

with b = (a + /3)/2 and c = d = {a — (3)/2. We recognize in ( [5.3| ) the partition function of 
the 8-vertex model on random dynamical lattices.H A configuration of the model is defined 
by a quadr-angulated surface with arrows on the edges of the graph, such that each of the 
vertices displays one of the eight allowed configurations, which are weighted with the 3 
constants 6, c, d. For a = (3 the t/(l)-breaking term X^ + X^^ vanishes and we recover 
the 6-vertex modelM 

We shall now show how to solve the model in the planar limit, i.e. compute the large 
A^ free energy. 



Or, more precisely, a two-parameter slice of it, since c and d are not independent. 

In this case note that the arrows "cross" each other just like knots in links, so that we 
manifestly recover our link model with the 2 orientations of the loops playing the same role as 
the 2 colors. 
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5.1. Character expansion 

All known matrix model solutions are (more or less implicitly) based on the fact that 
we can reduce the number of degrees of freedom from A^^ to N . Usually the N remaining 
degrees of freedom are the eigenvalues of the matrices. Unfortunately, from Eq. (|5.2| ) one 
cannot go directly to the eigenvalues of A and B: we do not know how to integrate over 
the relative angle between A and S.a Therefore, instead of working directly with (p. 21), we 



expand the troublesome part exp(A^tr(Ai?)^)) in characters of GL{N). Recalling that 
all class- functions can be expanded on the basis of characters, we write: 

iV|tr(AS)2 ^ 
e =5^C{^}X{/.}(^S) (5.4) 

{h} 

where X{h}{-^B) is the character taken at AB and {h} is the set of shifted highest weights 
hi = rui + N — i {rrii highest weights), hi > h2 > ■ ■ ■ > hn > 0, which parametrize the 
GL{N) analytic irreducible representation. The coefficients of the expansion c^/i} can be 
determined explicitly: 

.,., . (A.,/,)*V.^™|gl^i)Z!) (5.5) 

in terms of the set {/i^^*'"} and {h°'^'^} of even and odd hi. The advantage of characters 
is that they satisfy orthogonality relations, so that we can now integrate over the relative 
angle between A and B: 

I dnxw(-4nBfit).Ml(^)Ml(^ (5.6) 

JU{N) X{h}0-) 

where the dimension X{/i}(l) is up to an overall constant the Van der Monde determinant 
AN. 

Once Eqs. (|5.4| ) and ( |5.6|) are inserted into (|5.2|), we see that the integrand only 
depends on the eigenvalues of A and B: 



(^) (^, a\-^^ ^W 



ZABABi^^l3)=Yl 



A[hi 

{h} L 



N AT^N ' -" -.2 






(5.7) 



The key observation here is that we still have an action of order A^, but we have A 
highest weights hi and A eigenvalues A^; therefore a saddle point analysis applies again. 



4 



Only one integral of this type is known exactly, the Harish Chandra-Itzykson-Zuber integral 



[13,113], but it does not apply here. 
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5.2. Saddle point on Young tableaux 

The notion of a saddle point on Young tableaux first appeared in |jl^ in the context 
of the asymptotics of the Plancherel measure. It was rediscovered independently in the 
solution of large A^ 2D Yang-Mills [^ , and was used to deal with character expansions in 



21]. In the present calculation, the novelty is that we have to deal with a double saddle 



point equation on both eigenvalues and shifted highest weights (i.e. shape of the Young 
tableau) Q. 

The idea here is to find an appropriate scaling ansatz for the shape of the dominant 
Young tableau in the large N limit. We find that the highest weights hi scale as A^ 
(the Young tableaux become large both horizontally and vertically), so we can define a 
continuous density of rescaled hi/N: 

1 ^ 

i=l 

and the corresponding resolvent 

., Pih') 



H{h) = / dh 



h-h' 



We also have a density of eigenvalues p(A) and the resolvent uj{X). 

The saddle point equations now read (we use "slashed" functions defined by: ]^{h) : = 
^{H{h + iO) + H{h — iO)) and similarly for the other functions): 

-A + aA3+0(A)+^(A)/A = O Ag[-Ao,+Ao] 
l^(h)-^ = -log{h/P) he[h^,h2] 
with L{h) = logA^(/i). The new unknown functions h{X) and A(/i) appear when tak- 



h 



ing the logarithmic derivative of det^ ^[Aj^]; this type of functions was analyzed in P^ 



where it was shown that X{h) and h{X) are functional inverses of each other. There- 
fore, we have two saddle point equations which are connected by a functional inver- 
sion relation. This connection allows to solve them; skipping the details, one can 
show that one has a well-defined Riemann-Hilbert problem for the auxiliary function 
D{h) := 2L{h) — H{h) — 3 log h + \og{h — hi), whose solution can be expressed in terms of 
functions in an appropriate elliptic parametrization y{h): 

n/^,^ i h - hi log(/?/a) .,^ , „, e2{xo - y{h)) 

D{h) = log -^ y{h) + 2 log — 

—ah^ K ^2{xo + y{h)) 
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Fig. 11: Phase diagram of the ABAB model. The dashed hne is the a = (5 hne, 
the curves are equipotentials of the elUptic nome. 

5.3. Phase diagram and discussion 

In the same way that the one-matrix model displayed a singularity at Qc = 1/12, here 
the free energy and the various correlation functions have a line of singularities in the 
(a, /3) plane, which is shown on Fig. 11. 

We recognize at a = 1/12, /3 = the usual singularity of the one-matrix model. 
In fact, one can show that everywhere on the critical line except at the critical point 
Cic = (^c= 4^7 the critical behavior is the same as the one of the one-matrix model ("pure 
gravity" behavior). This implies the following asymptotics: at fixed slope s = P/a 7^ 1, if 
the free energy F{a,(3 = sa) = J2p fpi^)'^^ then 



fp{s) ~ 



const ac{s) ^ p 



-7/2 



However, the point etc = /3c = 4^, that is the 6- vertex model point, is very special: 
this is the point where the elliptic functions degenerate into trigonometric functions, which 
implies logarithmic corrections: 

/p(l) ^~°° const {ATiYp~^\ogp 

This is characteristic of a c = 1 conformal field theory coupled to gravity. 

5.4. Application to reduced alternating diagrams 

We should remember that the a = (3 line is of special interest to us, since it is the 
intersecting loops 0(2) model (solving a certain counting problem for alternating links) we 
started from. In order to carry out the program that we have applied to the one-matrix 



14 



model we should next address the two issues of primality/mininiality and of the flype 
equivalence. We shall only consider the first issue; a dicussion of the flype equivalence in 
the general 0{n) case will be made in the next section, and the corresponding calculation 
for n = 2 will appear in a future publication [^ . 



Again, we introduce an additional parameter t in the action: 

,^, f Ntr (-UA^ + B^) + f (A4 + B^) + UABf) 



and we impose that the 2-point function G2{t,g) — limjv^oo(;^trA^) — limiv^oo(;^trS^) 
satisfies 

G2{t{g),g) = l (5.10) 

Obvious scaling properties imply that G2{t, g) = |G2(1, d/t^): and the formula for 6*2(1, g) 
can be found in appendix A of [|11| in terms of complete elliptic integrals. This gives 



an equation for t[g), which can in principle be solved (at least to an arbitrary order in 
perturbation theory). 



In order to go further, we notice that at the singularity we must have gc/t{gc 



\2 _ 1 



47r' 



from [|ni we extract 6*2(1, l/(47r)) = |(4-7r), and therefore using j fU] ), t{gc) = f (4-7r), 
which finally yields 

gc = ^{n-4f (5.11) 

We conclude that the number fp of reduced alternating link diagrams with 2 colors and p 
crossings has the following asymptotics: 

where the number 1/gc = 6.91167... is slightly larger than the value 6.75 obtained for 
only one color. 



6. Further generalizations and prospects 

We have already written a fairly general model, the intersecting loops 0{n) model 
(Eq. (|3.2| )) which should contain in principle all information on the counting of alternating 
links and knots. We shall now show how this model is in fact not sufficient for our purposes. 

Indeed, one should remember that the 0{n) model given above counts alternating link 
diagrams, and not alternating links. For the latter, one should address the problem of the 
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flype equivalence. We have seen in the n = 1 case (section 4) that we needed to do a httle 
surgery on the four-point functions. One can convince oneself that what it amounts to, in 
more physical terms, is a finite renormalization which results in the appearance of quartic 
counterterms in the action. Generically these counterterms will have the most general 
form compatible with the symmetry. In our case, we find that there are two independent 
0(n)-symmetric tetravalent vertices, of the form MaMhMaMf, and MaMaM^Mh. This 
results in a generalized 0{n) model: 

f- ATtr (-1 ELi M! + f j:i,=iiMaM,r + I E:.=i M^^mA 

Z^''\n,g,h)= UdM^e ^ ^ 

•^ a=l 

(6.1) 
where h will be given as a function of g by appropriate combinatorial relations of the same 
form as those of section 4. 

At the moment, the solution of this general model is unknown. Let us note however 
that for /i = 0, this model is simply the usual (non-intersecting loops) 0{n) model which 
has been completely solved 0] . It is tempting to speculate that there is no phase transition 
in the ((/, h) plane as one moves away from the h = lineB Then, one can make predictions 
on universal quantities such as critical exponents of the model. For example, the number 
fp{n) of prime alternating links with n colors would have the asymptotics 

fp{n) ~ const(n) b{nY p~'^~^^'' n = -2cos(7rz/), < z/ < 1 

In particular the number fp of prime alternating knots would satisfy 

fp ^ const b{0)P p-^ 

One interesting question is to determine the non-universal constant 6(0). This, of course, 
requires to really solve the n = model (or more precisely to study the n ^ limit). An 
alternative option is to note that this model can be recast as a supersymmetric Osp{2n\2n) 
model, the simplest (n = 1) being a supersymmetrized version of the 0{2) model considered 
earlier; using bosonic and fermionic complex matrices X and ^, it can be written as: 

Z^^\g) = /"dXdXM^d^^ 

(6 2) 
A^tr (-XAt - ^^t + g{XX^X^X + ^^t^t^ + ^X^'H^X + At^X^t)) 



e 



This is certainly true in the n = 2 case, as the study of p3] shows, and the exact solution 



1 25] confirms. 
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Due to supersymmetry, this partition function is equal to 1, but non-supersymmetric 
correlation functions such as (-^tr(XX^)") are non-trivial and should contain the desired 
information. Whether this model is solvable or not is an open question. 
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